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The response of the QED vacuum in an asymptotically large electromagnetic field is studied. In this regime the vacuum energy is 
strongly influenced by the vacuum polarization effect. The possible interaction between the virtual electromagnetic radiation and 
a superstrong magnetic field suggests that a background of virtual photons is a source of magnetization to the whole vacuum. The 
corresponding contribution to the vacuum magnetization density is determined by considering the individual contribution of each 
vacuum polarization eigenmode in the Euler-Heisenberg Lagrangian. Additional issues concerning the transverse pressures are 
analyzed. We also study the case in which the vacuum is occupied by a superstrong electric field. It is discussed that, in addition to 
the electron-positron pairs, the vacuum could create photons with different propagation modes. The possible relation between the 
emission of photons and the birefringent character of the vacuum is shown as well. 
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1. Introduction 

Whilst there is some evidence that very large magnetic fields 
|B| » B c , B c = m 2 /e = 4.42 • 1O 13 G0 exist in stellar objects 
identified as neutron stars fl|,Q,|3|], its origin and evolution re- 
mains poorly understood Jj]. Some investigations in this area 
provide theoretical evidence that |B| might be generated due to 
gravitational and rotational effects, whereas other theories esti- 
mate that is self-consistent due to the Bose-Einstein condensa- 
tion of charged and neutral boson gases in a superstrong mag- 
netic field JUS SB- In this framework the nonlinear QED- 
vacuum possesses the properties of a paramagnetic medium and 
constitutes a source of magnetization, induced by the external 
magnetic field. Its properties are primarily determined by the 
vacuum energy of virtual electron-positron pairs. Because of 
this, a negative pressure transversal to the external field is gen- 
erated 1 9] in similarity with the Casimir effect between metallic 
plates llioll . Moreover, the vacuum occupied by the external 
field turns out to be an "exotic" scenario in which processes 
like photon splitting HH [H] and photon capture (ij, [3, Hi 
could take place. These two phenomena depend on the pho- 
ton dispersion relation which differs from the light cone, due to 
vacuum polarization effects 

iEE O, Ell As a result, the 

issue of light propagation in empty space, in the presence of an 
external magnetic field, is similar to the dispersion of light in 
an anisotropic "medium". 

The phenomenological aspects associated with this problem 
have been studied for a long time. In the meanwhile, other fea- 
tures of nonlinear electrodynamics in a superstrong magnetic 



1 Hereafter m and e are the electron mass and charge, respectively. 
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field have been studied such as the dimensional reduction of the 
Coulomb potential 1I20I.I21I.I22I.I23I1 and the possible existence of 
a photon anomalous magnetic moment 12411 . However, due to 
the vacuum polarization effect, virtual photons can carry a mag- 
netization as well. As a consequence, they might be a source 
of magnetism to the whole vacuum. Motivated by this idea, 
we address the question in which way the virtual electromag- 
netic radiation contributes to the vacuum magnetization and 
therefore to increase the external field strength. The magnetic 
properties of the vacuum have been studied in 
for weak (|B| «: B c ) and moderate fields (|B| ~ B c ) in one- 
loop approximation of the Euler-Heisenberg Lagrangian [ 28] 
which involves the contribution from virtual electron-positron 
pairs. The contribution of virtual photons, created and annihi- 
lated spontaneously in the vacuum and interacting with B by 
means of n^y, is contained within the two-loop term of the 
Euler-Heisenberg Lagrangian (see Fig. [TJ. The latter was com- 
puted many years ago by Ritus I 29i 30n and has been recalcu- 
lated by several authors as well [EnBlIISI. In all these 
works, however, it is really cumbersome to discern the indi- 
vidual contributions given by each virtual photon propagation 
mode to the Euler-Heisenberg Lagrangian which should allow 
to determine the magnetism and pressure associated with each 
form of virtual mode. In this Letter we analyze these contribu- 
tions separately for very large magnetic fields (|B| » B c ) since 
these allow to establish relations between the birefringence of 
the vacuum lfl9ll35l.l36ll and the global properties of it. 

Besides the strongly magnetized vacuum, there is another 
interesting external field configuration which deserves to be an- 
alyzed: a supercritical electric field |E| » E c with E c = m 2 /e = 
1.3 • 10 I6 V/cm. In this asymptotic region the Euler-Heisenberg 
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Figure 1 : Two-loop expansion of the Euler-Heisenberg Lagrangian. The double 
lines represent the electron-positron Green's functions, whereas the wavy line 
refers to the photon. Here X* ' is the free Maxwell Lagrangian, X*" represents 
the one-loop which gives the contribution of the virtual free electron-positron 
pairs created and annihilated spontaneously in vacuum and interacting with the 
external field. The radiative corrections (involved in X' 2 ') emerge from two- 
loop due to exchange of the virtual photons. 



Lagrangian acquires an imaginary term which characterizes the 
instability of the vacuum. This phenomenon is closely related 
to the production of observable particles from the own vacuum. 
Certainly, the creation of electron-positron pairs -the so-called 
Schwinger mechanism- turns out to be the most remarkable ef- 
fect predicted through this procedure J28, 37, 3^1 . However, 
the imaginary part of this effective Lagrangian is just a mea- 
sure of the vacuum decay and does neither give the actual rate 
of production of particles nor the accessible decay channels 
1 39]. Thereby not only the creation of electron-positron pairs 
is a plausible effect but also the emission of observable pho- 
tons 1 29[ 3(J] • The latter phenomenon was analyzed by Gitman, 
Fradkin and Shvartsman IdOIdII 15211 . Their results showed that 
the total probability of photon emission from the vacuum, ac- 
companied by the creation of an arbitrary number of electron- 
positron pairs, is connected to the decay probability of the vac- 
uum and thus to the imaginary part arising from the two-loop 
term of the Euler-Heisenberg Lagrangian. In this context the 
corresponding decomposition in terms of the vacuum polariza- 
tion modes is particularly illuminating because it reveals that 
only two of them contribute to the vacuum instability. It seems, 
therefore, that the vacuum could create photons with different 
propagation modes, an effect closely related to its own birefrin- 
gence. 



2. Preliminary remarks 

In a magnetized vacuum the spatial symmetry is explic- 
itly broken by the external field B. In this context, there is a 
vectorial basis b? Ir20[ 4oL 41] which characterizes the vac- 
uum symmetry properties and fulfills both the orthogonality 

condition: b^b^ = S'j\b®) and the completeness relation: 

d>* v - = 2i=i ^bvVfb®) 2 - Explicitly, the basis vectors 
read b^ = k 2 .^ 2 k A - k^k& 2 k), bf = &^k x , bf = ^ A k A 



and b^ 4) = k^. These expressions involve the external field ten- 
sor J\, v and its dual # /JV 
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vacuum polarization tensor is diagonal i.e. 

4 l(0u© 

U (b (!) ) 2 



In this basis, the 



(1) 



whereas the dressed photon Green function can be expressed as 
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k 2 k 2 
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Here the represent the TL^ -eigenvalues and £ is the gauge 
parameter. This diagonal decomposition of n^y defines the en- 
ergy spectrum of the electromagnetic field which differs from 
the isotropic vacuum (B = 0). 

Owing to the transversality property (kTLj V = 0), the eigen- 
value corresponding to the fourth eigenvector vanishes identi- 
cally (xr (4) = 0). Furthermore, not all the remaining eigenmodes 
are physical. In general, this depend on the direction of wave 
propagation. To show this we consider b? (&) as the electromag- 
netic four vector describing a photon. The corresponding elec- 
tric and magnetic fields of each mode are e w = i(kb® - o»®b®), 
ft w = —ik x b w . It follows that the mode i = 3 is a wave polar- 
ized in the transverse plane to k whose electric e (3) ~ k ± x «y 
and magnetic Z> (3) ~ n\\k\ - kjt\\ fields are orthogonal to B as 
well. Here the vectors k ± and k\\ are the components of k across 
and along B with n\\ = B/|B|. For a pure longitudinal propaga- 
tion to the external field k ± = 0, the mode is a longitudinal 
and non-physical electric wave e <2) ~ n\\. On the other hand, 
b^, is transverse since the associated electric field is e (1) ~ k ± 
whereas the magnetic one b m ~ k ± X k\\. As a consequence, 
both b^ 1 ' and bf represent physical waves which may be com- 
bined to form a circularly polarized transversal wave. In this 
case both modes propagate along B with a dispersion law inde- 
pendents of the magnetic field strength B17l ll9. 41]. 

Now, if the photon propagation involves a nonvanishing trans- 
versal momentum component k± + 0, we are allowed to per- 
form the analysis in a Lorentz frame that does not change the 
value k ± , but gives k\\ = and does not introduce an external 
electric field. In this Lorentz frame, the first eigenmode b^ 1 ' 
becomes purely electric longitudinal and a non physical mode 
whereas b^ is transverse. Hence, for a photon whose three- 
momentum is directed at any nonzero angle with the external 
magnetic field, the two orthogonal polarization states b^, 2) and 
b^ 3) propagate. In this framework the analytical structures of 
the corresponding eigenvalues 2*2,3 are different. As a matter 
of fact, the vacuum behaves like a birefringent medium with 
refraction indices II 19L l4lh 
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Here 0)2,3 are the corresponding solution of the dispersion equa- 
tions k 2 = arising from the poles of 

Considering these aspects, we analyze the Euler-Heisenberg 
Lagrangian 

(4) 



Xeh - X-r + X.R + ■ 



where jj^ = -|B 2 is the free renormalized Maxwell Lagrangian, 
whereas jj^ denotes the one loop regularized contribution of 
virtual electron-positron pairs created and annihilated sponta- 
neously in vacuum and interacting with B 1 28 ] . In asymptoti- 
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Figure 2: Diagrammatic decomposition of X <2 ' in terms of the vacuum polar- 
ization eigenmodes. 



cally large magnetic fields it reads 
m 4 fa 2 



24n 2 



In 
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Here fa = |B|/B C , ln(y) = 0.577 . . . denotes the Euler constant 
whereas 6n~ 2 £'(2) = -0.5699610 . . . and g(x) is the Riemann 
zeta-function. 

The contribution of a virtual photon interacting with ex- 
ternal field by means of the vacuum polarization tensor is ex- 
pressed as S (HQ 

£(2> = 5 / T^n^k)^), (6) 

where v (k) denotes the bare photon Green function which is 
obtained by neglecting x,- in the denominator of f^ v . Using this 
expression we find 
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wherecffc = d k/(2n) '. To obtain the above expression we have 
inserted the diagonal decomposition of Tl^ into Eq. © and 
used the orthogonality condition. Certainly, the quantity Xjk~ 2 
represents the interaction energy density of a virtual photon 
propagation mode with B. Because the analytical properties of 
the Hi differ from each other the contribution of each mode will 
be different. Therefore, the original two-loop graph can be de- 
composed into three diagrams (see Fig. [2]i whose properties are 
connected to the birefringence property of the vacuum. There- 
fore, each individual term £?' = j Jdk>Ci/k 2 will be studied 
separately. As we shall see very shortly, the decomposition of 
£ 2) is also valid in an external electric field in which case some 
evidence on the possible decay of the vacuum into observable 
photons emerges. 

Obviously, the renormalized two-loop term decomposes as 



well: 



£r > = Hl=i -£;r aQ d tne leading behavior of £^ 



asymptotically large magnetic field is given by (for details see 
Ref. 
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These expressions are accurate up to terms that decrease with fa 
like ~ b ln(fa) and faster. Here the numerical constants are 



M = 1.25, N 2 = 0.71, N 3 = 0.29, 
JK = 4.21 and S = 0.69. 
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Moreover, the asymptotic behavior of the full two-loop term is 
am 4 fa 2 
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[yn) 



+ 2.4 



(10) 



which coincides with the results reported in references |2§l 31]. 



3. Properties of a highly magnetized vacuum 

In the presence of an external magnetic field, the zero point 
vacuum energy <? vac is modified by the interaction between B 
and the virtual QED-particles. The latter is determined by the 
effective potential coming from the quantum-corrections to the 
Maxwell Lagrangian which is also contained within the finite 
temperature formalism. According to Eq. © it is expressed as 



<?vac = - -£ R 1)- 2f=i £;r + - ■ ■ Consequently the vacuum acquires 
a non trivial magnetization ^# V ac = -d^ vac /9\B\ induced by the 
external magnetic field. In what follows we will write ^# vac = 
■/#vac + ^vac + • • ■ in correspondence with the loop-term £^. 
In this sense, the one loop contribution at very large magnetic 
field b » 1 can be computed by means of Eq. (|5]l and gives 
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Incidentally, the above asymptotic behavior is also manifest in 
the corresponding magnetization derived from the QCD vac- 
uum in a magnetic field |B| » AQ CD /e. For more details we 



refer the reader to Ref. J42 ] . 

The two-loop correction is given by ,/# (2) = Ya=i ■ 
where ./^f 2) = d£^/d\B\ is the contribution corresponding to a 
photon propagation mode. Making use of Eqs. © we find 
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(12) 
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According to these results, in a superstrong magnetic field 
approximation, jtf^ < and ^# 2 (2 ' < ^ behave diamagnet- 
ically whereas j^T* > is purely paramagnetic. Moreover, 
while ^#^ 2) depends linearly on fa, the contributions of the sec- 
ond and third propagation mode depend logarithmically on the 
external field. We find, in particular, that for magnetic fields 
B ~ 10 18 G, the magnetization generated by the first and second 
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polarization mode reaches values the order ~ -10 12 erg/(cm 3 G) 
and ~ - 10 13 erg/ (cm 3 G), respectively. In the same context ^# (3) 
+ 10 13 erg/(cm 3 G). Note that the leading behavior of the com- 
plete two-loop contribution is 



am b 
167r 3 B c 



ln| — | + 2.9 

yn 



>0. 



(15) 



which shows a dominance of the third mode. Indeed, for B ~ 
10 18 G, one finds Jt {1) ~ +10 12 erg/(cm 3 G). 

As it was expected ^# (1) /^ (2) ~ ce~ x . This ratio is also 
manifested between the corresponding magnetic susceptibili- 
ties (JT® = d^# ( ''V<9|B|). Note that 
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+ 1.86 



> 0, 
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(16) 
(17) 



For magnetic fields b ~ 10 5 corresponding to |B| ~ 10 18 G, the 
magnetic susceptibility reaches values of the order of ~ 
10~ 4 erg/(cm 3 G 2 ) which exceeds the values of many laboratory 
materials, for example Aluminum (^ai = 2.2T0~ 5 erg/(cm 3 G 2 )). 

Some additional comments are in order. First of all, even 
though the previous decomposition of ./# (2) is not really ob- 
servable, it turns out to be a transparent framework which illus- 
trates, in a phenomenological way, the magnetic property gen- 
erated by each photon propagation mode and thus, a connection 
with the birefringent property of the vacuum. The decomposi- 
tion of ..# (2) in terms of the photon modes allows, in addition, 
to establish similarities and differences with the magnetization 
carried by observable photon^ Indeed, similar to the latter the 
virtual radiation carries a magnetization which depends on the 
polarization vector. However, while an observable second mode 
has a paramagnetic response [4-5], our result points out that the 
corresponding virtual polarization generates a purely diamag- 
netic magnetization. Besides, we have seen that the first propa- 
gation mode contributes to the magnetization of the vacuum. 
This is not expected for an observable mode-1 photon since 
its dispersion law is independent on the external field strength 
112, 19, 4l|] and therefore does not carry a magnetization. 

Because of the anisotropy generated by B a magnetized vac- 
uum exerts two different pressures Q5|, |a |t1 W ■ One of them is 
positive {&\\ = -^vac) and along B, whereas the remaining is 
transverse to the external field direction (3 g ± = - ^|B|). 
For b ~ 1 the latter acquires negative values. At very large 
magnetic fields (b » 1) the one-loop approximation of & ± can 
be computed by making use of Eq. (0 and Eq. ( fTTY In fact 
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-The magnetic response of an observable photon was studied in 1 45] in two 
different regimes of the vacuum polarization tensor. On the one hand for low 
energies in weak fields (|B| « B c ) and on the other hand (originally studied 
in Ref. 1 44]) near the first pair creation threshold and for a moderate fields 

(|B| ~ B c ). 



Therefore, at asymptotically large values of the external field, 
the interaction between B and the virtual electron positron pairs 
generates a negative pressure which would tend to shrink in- 
serted matter in the plane transverse to B. 

Again, the two-loop contribution can be written as the sum 
of the corresponding terms due to the vacuum polarization modes 
Hii ^ff- According to Eqs. ® and Eqs. (O they 
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with the complete two-loop term given by 
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(20) 



(21) 



(22) 



For b ~ 10 5 , corresponding to magnetic fields B ~ 10 18 G, 
the transverse pressure generated by the first and second po- 
larization mode are positive and reaches values of the order 
~ 10 30 dyn/cm 2 and ~ 10 31 dyn/cm 2 , respectively. In contrast, 
the contribution given by the third mode is negative with ~ 
-10 31 dyn/cm 2 . The combined result is negative and achieves 
values of the order ^f 1 10 31 dyn/cm 2 . This fact strength- 
ens the analogy of the considered problem with the Casimir 
effect in which the pressure transversal to the parallel plates is 
also negative and dominated by the virtual electromagnetic ra- 
diation 04311 • 



4. The vacuum instability in a supercritical 
electric field 

Ultra-high electric fields |E| » E c = m 2 /e = 1 .3 • 10 I6 V/cm 
have been predicted to exist at the surface of strange stars (0 
47, 48, 49| . In this electric field regime, the asymptotic behav- 
ior of £^J? is obtained from Eqs. © by means of the duality 
transformation b — > -it with e = |E|/E e . As a consequence we 
can write 

XEH = ^[XEH]+Im[£ EH ] (23) 



where !R[Xeh] constributes to the dispersive effects. Because 
of the imaginary part the vacuum becomes unstable and cre- 
ation of particles could take place. The probability associated 
with the vacuum decay is P = 1 - |<0 out |0i n >| 2 with (0 out |0i n ) = 
giVTX E H ; w here VT is the volume element in 3 + 1 dimensions. 
With this in mind and by considering Eq. ( 1231 one has 
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y _ e -2Im[X EH ]VT 



(24) 
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As we have already mentioned in the introduction, the emis- 
sion of observable photons from the vacuum is also a plausible 
effect. The probability *p of photon emission from the vacuum 
accompanied by the creation of an arbitrary number of electron- 
positron pairs, can be determined by using the unitarity condi- 
tion for the dispersion matrix 5=1+ iSF which leads to a 
relation of the optical-theorem type lf50L 51 , 521 



term is positive 



Yj Kout|5Hin>| 2 = 2Im<in|^|in), 



(25) 



where 3" is the sum of all Feynman graphs. However, the elec- 
tron causal Green function for the matrix elements (in| . . . |in) 
is quite different from the standard propagator in the Schwinger 
proper- time representation of the out-in matrix elements. We 
consider the relation above with the in-state having no pho- 
tons i.e. |in) = |0j n ). In addition, we will confine ourselves 
in the right- and left-hand sides of Eq. ( f25l l to the second or- 
der radiational interaction. In this approximation, one can take 
the operator 3 on the left-hand side to first order: S?^ - 
- J d 4 x j^(x)a ll (x). Here f = | [0 , (jc)y A ', t//(x)] is the current and 
a^ix) denotes the radiation field. Therefore, only one-photon 
states contribute to the sum over the out-states (out| = (Oi n |ct,- 
b 9] . . . b qi d Pi . . . d Pr Whilst c^, is understood as the annihilation 
operator of a mode-/ photon, b and d are interpreted as the an- 
nihilation operators of electrons and positrons, respectively. On 
the right-hand side of Eq. d25l l 3 has to be taken in second 
order 3 {1) = -| J d 4 xd 4 x'T [f(x)f(x , )a ft (x)a v (x')\ , where T 
represents the time ordering operator. Considering the normal 
mode expansions of the dynamical fields and by using the Wick 
theorem one obtains 



(in|^ (2) |in> = £ { 



(2) 



(26) 



where is determined by Eq. (O with the causal Feynman 
propagators replaced by the corresponding electron Green func- 
tion appearing in the mean values <in| . . . |in) of Eq. (l25l l. 

Now the decay probability P is connected to the total prob- 
ability of photon emission from the vacuum, accompanied by 
the creation of an arbitrary number of electron-positron pairs 
P = *p + . . . Obviously, the corresponding expansion upto 
second order in a involves the imaginary part coming from 
the two-loop contribution of the Euler-Heisenberg Lagrangian. 
The latter can be written as the sum of the corresponding terms 
due to the individual vacuum polarization modes Im [-£^j = 
^ =1 Im[4^]with 



(27) 



Note that, because of the dominance of the third propagation 
mode, the complete imaginary part arising from the two-loop 



An 



(28) 



3 For details, see the cited references of Fradkin, Gitman, Shvartsman. 



where Im [-£^j - is the imaginary part coming from the 
one-loop contribution at very large electric field e » 1 . 

There is some interesting aspects in Eq. ( 127b which de- 
serves some comments. First of all, it reveals that the first prop- 
agation mode does not contribute to the imaginary part of Xeh- 
This suppression provides an evidence that the vacuum does 
not generate electromagnetic waves propagating along the ex- 
ternal field E. On the other hand, only the second and third 
mode contribute to the vacuum instability. This is a signal that 
the vacuum could create the corresponding propagation modes. 
However, the birefringence character of the vacuum imposes 
that there is not a single observable mode-2 photon without the 
existence of a corresponding mode-3 photon. Therefore, in ad- 
dition to the electron-positron pairs, the vacuum could create 
photons with different propagation modes. Note, however, that 
the described photon emission has been predicted within the 
framework of equilibrium quantum field theory, even though it 
is a far-from-equilibrium, time-dependent phenomenon. Only 
further studies can tell us how far this process can be stretched 
because a realistic treatment of this issue requires a time evolu- 
tion analysis of the photon number distribution functions, simi- 
lar to that developed by Hebenstreit et al. for electron-positron 
pairs within a quantum kinetic approach ll53ll . 

5. Summary and Outlook 

In summary, we have examined the magnetization of the 
QED vacuum in the presence of a constant magnetic field in 
the strong field regime, |B| » B c . We have seen that the virtual 
electromagnetic radiation is a source of magnetization to the 
whole vacuum. In a superstrong magnetic field approximation, 
the two-loop contribution of the magnetization density corre- 
sponding to the second and third propagation mode depends 
nonlinearly on the external magnetic field and their behavior is 
diamagnetic and paramagnetic, respectively. On the other hand, 
the contribution coming from the first mode is diamagnetic and 
depends linear on B. We have seen that for very large magnetic 
field the contribution of the third mode strongly dominates the 
analyzed quantities. In this regime the latter tends to shrink in- 
serted matter by exerting a negative transverse pressures to the 
external field. On the contrary those contributions coming from 
the first and second virtual mode are positive and tend to expand 
the matter. 

In the last section of this Letter we showed that only two 
photon propagation modes contribute to the instability of the 
vacuum in an strong electric field. This instability is associated 
with the emission of photon whose propagation modes differ 
each other. A plausible connection between this mechanism 
and the birefringent character of the vacuum occupied by a su- 
percritical electric field was discussed and the suppression of 
vacuum decay into pair of modes propagating along an external 
electric field was analyzed as well. 
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